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We show that puncture data for quasicircular binary black hole orbits allow a special gauge choice 
that realizes some of the necessary conditions for the existence of an approximate helical Killing 
vector field. Introducing free parameters for the lapse at the punctures we can satisfy the condition 
that the Komar and ADM mass agree at spatial infinity. Several other conditions for an approximate 
Killing vector are then automatically satisfied, and the 3-metric evolves on a timescale smaller than 
the orbital timescale. The time derivative of the extrinsic curvature however remains significant. 
Nevertheless, quasicircular puncture data are not as far from possessing a helical Killing vector as 
one might have expected. 
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I. INTRODUCTION 

Binary black hole (BH) mergers are among the most 
promising sources for the ground-based interferomet- 
ric gravitational wave detectors GEO600, LIGO and 
TAMA U , which have started to collect data. Predicting 
the gravitational waves from the dynamic and non- linear 
regime of BH mergers requires numerical simulation, and 
a crucial issue is to find astrophysically realistic initial 
data. This is a non-trivial task which involves specifica- 
tion of some appropriate free data and the solution of the 
constraint equations of general relativity for the remain- 
ing data. By now several methods to produce initial data 
for binary BH systems exist (see for a recent review). 
However, there are still open issues concerning the con- 
struction of astrophysically realistic initial data, which in 
particular are data that represent two BHs on almost cir- 
cular orbits during their inspiral, and furthermore data 
that represent a quantifiable approximation to two BHs 
near the transition from inspiral to plunge and merger. 

During the inspiral, we expect the two BHs to be 
in quasicircular orbits around each other with a radius 
which shrinks on a timescale much larger than the or- 
bital timescale. This means that the initial data should 
have an approximate helical Killing vector £. In addi- 
tion one would like to have the initial data in coordinates 
such that this approximate symmetry is manifest, i.e. the 
time evolution vector should lie along £, so that the time 
derivatives of the evolved quantities are minimized. Then 
we would get 

d t9ij w d t <t> « d t K a d t A l3 a 0, (1) 

where we have decomposed the metric into a conformal 
factor <f> and a conformal metric gij, and split the ex- 
trinsic curvature into its trace K and a tracefree piece 
Aij. All the terms in Eq. are of the order of some 
error quantity which measures the failure of an exact he- 
lical Killing vector to exist. This error has a finite value 
near the "innermost stable circular orbit" and tends to 



zero for increasing separation of the BHs. In principle 
one can improve the approximation by replacing the cir- 
cular orbits by inspiraling orbits defined by a small but 
nonzero value of all the time derivatives in Eq. (yi based 
on post-Newtonian approximations, see for example 

One approach to construct quasicircular initial data 
is to use the conformal thin sandwich (CTS) decompo- 
sition 0, ||. For CTS data the time derivative of the 
conformal metric is free data that can be set to zero. If 
in addition a maximal slicing lapse is used, one obtains 
d t gij — d t K = 0. Note, however, that d t (j) and d t Aij are 
in general non-zero, since CTS by itself does not contain 
conditions concerning helical Killing vectors. In the case 
of the Meudon data [E0] (but see also H@), CTS data 
is constructed under the assumptions of a two-sheeted 
topology with isometry, a conformally flat 3-metric, and 
a vanishing K . In addition, in order to construct qua- 
sicircular orbits, boundary conditions on the lapse and 
shift are imposed such that the Arnowitt-Deser-Misncr 
(ADM) and Komar mass (computed from the lapse) are 
equal, which is a necessary condition for the existence of 
a helical Killing vector. The expectation is that these 
conditions yield 

d t9l3 = dtK = 0, d t <t> « d t A l3 a 0. (2) 

It is not clear a priori how well dt<j> ~ and dtAij a 
are satisfied since only one necessary condition for the 
existence of the helical Killing vector has been enforced, 
namely the equality of the ADM and Komar mass. It 
would be interesting to check how much dt<f> and dtAij 
deviate from zero. Compared to Eq. l[T)l. there are now 
additional errors in Eq. (0, which are related to the 
assumption of conformal flatness, and to the details of 
the construction. In particular the exact imposition of 
dtgij — and the method used for determining the or- 
bital angular velocity of quasicircular orbits may intro- 
duce extra errors. 

Another approach to construct binary BH initial data 
are puncture data [ljj, which are calculated using the 
conformal transverse traceless (CTT) decomposition 0, 
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In this decomposition dtgij = cannot be achieved 
in general, while d t K = can always be imposed if a 
maximal slicing lapse is used. There still remains the 
gauge freedom to choose a shift that makes some of the 
time independence of a helical approximate Killing vec- 
tor manifest. So far approximate helical Killing vectors 
have not explicitly been used in conjunction with punc- 
ture data. Instead, in order to obtain quasicircular orbits 
the effective potential method E2 was use d m the past. 
In this method quasicircular orbits are identified with the 
cxtrcma in the binding energy. For thin spherical shells 
of point particles it can be shown 0] that both the effec- 
tive potential method as well as equality of Komar mass 
and ADM mass lead to the same identification of circu- 
lar orbits. In general, however, it is not clear whether 
the effective potential method leads to an approximate 
Killing vector in some sense. 

The angular velocity Q at the innermost stable circu- 
lar orbit (ISCO) computed from CTT puncture data with 
the effective potential method is almost twice as large as 
for Meudon data. Since the fl found in the Meudon data 
is close to post-Newtonian predictions, the CTT method 
may be less reliable for ISCO data. It is important to 
note, however, that if quasicircular CTT data is com- 
puted for the separation predicted by the post-Newtonian 
and the Meudon CTS methods, then the discrepancy in 
fl amounts to only a few percent |l4j | . This is an indica- 
tion that it is not the CTT method by itself but rather 
the particular definition of the ISCO that is problem- 
atic. The remaining discrepancy may have to do with the 
different conformal decompositions used, but could also 
come from different boundary conditions, or the differ- 
ent ways in which circular orbits are found. The key dif- 
ference between quasicircular data in the CTS and CTT 
methods is probably the construction of the extrinsic cur- 
vature. For CTS it is chosen such that dtgij = 0, while 
for CTT the Bowen-York extrinsic curvature is used with 
the hope that the orbital parameters can be adjusted such 
that reasonable approximations to circular orbits are ob- 
tained. An analytical comparison of CTT data and CTS 
data can be found in [l5|. 

In this paper we focus on CTT puncture data, for 
which we want to mention three reasons. First, the 
discrepancy between certain CTS and CTT approaches 
is not fully understood, but CTT appears to be more 
problematic. Second, in principle the puncture construc- 
tion is significantly simpler than current implementations 
of CTS with excision. And third, to date all gravita- 
tional wave forms obtained numerically for binary BH 
inspirals are based on puncture initial data, see e.g. 

El IB 13 El 13. 

In the following, we investigate the quasicircularity of 
CTT puncture data numerically by asking whether such 
data possesses an approximate Killing vector. In our 
analysis we do not make changes to the standard CTT 
puncture data construction, but we try to find a lapse 
and shift such that Eq. Q holds approximately. Con- 
cretely, we first compute a particular puncture data set 



for a quasicircular orbit as defined by the effective poten- 
tial method. We then compute a maximal slicing lapse 
with novel boundary conditions that are derived from 
necessary conditions for the existence of a Killing vector, 
namely the equality of the Komar and ADM masses at 
infinity or the punctures. Given the puncture data and 
the lapse, we calculate a shift which minimizes dtgij in 
the sense that a certain divergence of the metric does not 
evolve in time, obtaining 

d t d j9 11 = 8 t K = 0, d tgi j » d t <f> w 8 t A l3 « 0. (3) 

In other words, we use the four degrees of freedom 
in lapse and shift to set the time derivatives of four 
metric and curvature quantities to zero. Compared to 
Eqs. and J5J, in addition to the assumptions about 
the Killing vector and conformal flatness there will be 
non- vanishing terms due to the choice of Bowen-York ex- 
trinsic curvature and the effective potential method. We 
therefore check numerically how small the time deriva- 
tives are that are supposed to vanish approximately in 
Eq. ©, in particular we use the Baumgarte-Shapiro- 
Shibata-Nakamura (BSSN) system [2ll |22| of evolution 
equations to quantify the magnitude of the time deriva- 
tives. 

The paper is organized as follows. In Sec. [n] we de- 
scribe puncture data and briefly explain the methods we 
use to construct them numerically. Sec. lIIII describes how 
the helical Killing vector assumption can be used to con- 
struct a lapse and shift for puncture data and presents 
numerical results. In Sec. IIVI we conclude with a discus- 
sion. 

Notation. In this paper we use the standard 3 + 1 
decomposition of Einstein's equations, in which the 4- 
metric is written as 

ds 2 = -o?dt 2 + gij (dx i + I3 l dt) {dx 1 + I3 j dt) . (4) 

Here gij is the intrinsic 3-metric on a t = const hyper- 
surface and a and f3 z denote lapse and shift. Spatial 
components of tensors are denoted by Latin indices from 
the middle of the alphabet (e.g. i, j, ...). Latin letters 
from the beginning of alphabet are used for spacetime 
tensor indices. The extrinsic curvature is defined by 

K a b — ~-^£ n g a b: (5) 

where n a is the unit normal to the t = const hypersur- 
face. This implies that the time evolution vector is given 
by 




The 3-metric is decomposed into a conformal factor <j> 
and a conformal 3-metric g t j , such that 

9ij=<^9ij- (?) 
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The extrinsic curvature Kij is split into its trace K and 
its tracefree part A 1 ^ by writing it as 



(8) 



As a rule, physical quantities have an overbar while con- 
formal quantities do not, but lapse and shift are denoted 
by a and (3 l and are not rescaled. 



II. PUNCTURE DATA AND THE CTT 
DECOMPOSITION 

Puncture data are computed using the conformal 
transverse traceless (CTT) decomposition. In the CTT 
decomposition A VJ is decomposed into 



A 1 - 



+ LW %3 



(9) 



where A l ^ T is an arbitrary transverse traceless piece (i.e. 
S/jA l J, T = 0), and 



lw 1] = vr + V j W i - \g ij V k W k (10) 

is the longitudinal piece. Then the Hamiltonian and mo- 
mentum constraints become 



v- 



-r 7 {A^ T + LW^){A^ T + LW a )g lkgjl = 0(11) 



and 



V,LW lj 



a6v7» 



V l K = 0, 



(12) 



which are elliptic equations for <f> and W l . 

In the case of puncture data, the above CTS equations 
are further simplified by assuming that the 3-metric is 
conformally flat, i.e. 



In addition we set 



9ij = 4> 4 $i 



K = 



(13) 



(14) 



so that the slice is maximal. We also assume that the 
transverse piece vanishes, 

A% T = 0. (15) 

Analytic solutions to Eq. (|12J) are then given by 

Wi =E[-4^( 7P A + s Wl)l (16) 



A=l 



and hence 
K ij = A^ = cf)- w LW lj 

2 



A=l 



(17) 



fulfills the momentum constraint. The notation here is 
as follows: The coordinate locations of the two BHs are 
denoted by (x%,yx, z\) and (x 2 ,y2, z-i) and we have intro- 
duced 



r A = V(x - x A ) 2 + (y- va) 2 + (z- z a ) 2 . (18) 



and 



s\ = (x - x A , y-y A ,z- z A )/r A , 



(19) 



where the subscript A labels the BHs. The parameters 
P\ denote the momentum of each BH, and the ADM 
momentum at infinity is P^ + P^- On the other hand, the 
ADM momentum at the punctures, which represent the 
inner asymptotically flat ends of the hypersurface, van- 
ishes by construction, i.e. the black holes 'do not move' 
when viewed from the inner ends which is a useful feature 
when looking for an approximate Killing vector. Choos- 



ing PI 



we obtain data in coordinates where 



the net linear momentum vanishes at all three spatial 
infinities. 

In order to find a conformal factor which fulfills the 
Hamiltonian constraint we make the ansatz 



1 



mi 
2^ 



m 2 
2^ 



u, 



(20) 



which together with Eq. I|ll|) leads to 



V 2 u + \^- 7 LW^LW kl g lk g 3l = 0, (21) 



and we impose 



lim u = 

r — >oo 



(22) 



as boundary condition on u. The equation for u has 
to be solved numerically. The solution u is finite and 
at least twice differentiable. It depends on the bare 
masses m Al the momenta P A and the separation D = 
y/ (xi - x 2 ) 2 + (yi - V2) 2 + (zi - z 2 ) 2 . 

In order to obtain quasicircular orbits, we use the pa- 
rameters of Tab. IB The column denoted ISCO corre- 
sponds to the innermost stable circular orbit (ISCO) as 
determined by Baumgarte |23| with numerical values as 
in [24||. The Pre-ISCO column is a stable circular orbit 
found by Cook ^3 transcribed into puncture data [l4j . 
Note that for ease of comparison we use the same values 
for D, \P A \, and m A as in 0,13' which were normal- 
ized such that the ADM mass at infinity should be unity. 
However, we find that M^ DM deviates from unity by 
roughly 1%. We thus estimate these parameters to be 
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TABLE I: Parameters used in order to obtain quasicircular 
orbits within the effective potential method, with Q e // denot- 
ing the inferred angular velocity. The column denoted ISCO 
corresponds to the innermost stable circular orbit (ISCO) as 
determined by Baumgarte, and the Pre-ISCO column corre- 
sponds to a stable circular orbit found by Cook. The Pa = 
column is Brill-Lindquist data, i.e. two punctures with no mo- 
mentum. 



Parameter set 


ISCO 


Pre-ISCO 


Pa = 


D 


2.303 


3.698 


3.698 


\Pa\ 


0.3350 


0.2148 





tua 


0.4500 


0.4775 


0.5 


M ADM 
lvJ -oo 


1.003 


1.013 


1 


j ADM 


0.7715 


0.7944 







0.176 


0.102 






accurate to within about 1%, meaning that according to 
the effective potential method our orbits are not exactly 
circular. We were not able to obtain more accurate pa- 
rameters from the literature. 

Let us briefly discuss our numerical method. We solve 
the elliptic equations numerically using second order fi- 
nite differencing and a multigrid elliptic solver. The 
code is a stand-alone version of BAM_Elliptic [25|. All 
grids have uniform resolution. Unless explicitly speci- 
fied, all the results shown below are obtained with the 
outer boundary located at 12 and a finest resolution of 
h = 0.0625. As outer boundary conditions we use Robin 
conditions for all scalars, i.e. we assume that u ex v oc 1/r, 
where r is the distance to the center of mass. Here v de- 
termines the maximal slicing lapse (see Eq. (|31H below). 
As boundary condition for the shift (3q (computed below) 
we use /3q oc 1 jr 2 , which is a simplifying assumption that 
works reasonably well in practice. Since we are dealing 
with punctures, no inner boundaries are present. For the 
numerical work in this paper we consider non-spinning 
equal mass binaries with their center of mass at rest at 
the origin. The two BHs are positioned on the y-axis 
and their momenta point in the positive and negative x- 
directions, resulting in an angular momentum along the 
z-direction. After solving for u and v we observe second 
order convergence to zero in the Hamiltonian constraint 
and in dtK, as expected. We also find second order con- 
vergence in the shift f3 l after solving Eq. (|55|) below. The 
ADM and Komar integrals (denned in Sec. IIIII below) at 
infinity are computed with volume integrals over u and v 
covering the numerical domain plus a correction for the 
missing volume such that the overall error is expected to 
fall off like one over the distance to the outer boundary 
squared. At the punctures, ADM masses and Komar in- 
tegrals are obtained by fourth order interpolation of u 
and v onto the location of the punctures. From Fig. 2] we 
see that for h = 0.0625 the error in our masses at infinity 
is about 0.0005%, while the ADM and Komar integrals 
at the punctures have errors of 0.015% and 0.030% re- 
spectively. 
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FIG. 1: Shown (from top to bottom) are the ADM mass 
at infinity, the ADM mass at one of the punctures and the 
Komar mass at one of the punctures, each at three resolutions 
h for different locations of the outer boundary, in the case of 
the Pre-ISCO parameter set. The error in the ADM mass 
at infinity at resolution h — 0.0625 is about 0.0005%, while 
the ADM and Komar masses at the punctures have errors of 
0.015% and 0.030%, 



III. HELICAL KILLING VECTORS 

After solving for u, the puncture data are completely 
determined while lapse and shift are still arbitrary. We 
want to choose lapse and shift such that there is as little 
dynamical evolution as possible. This means we want a 
lapse and shift which satisfy 



d t g i: j = d t 4> = d t K = d t A^ = 



(23) 



as well as possible. If Eq. IL'MI) really did hold, it would 
imply the existence of a Killing vector 



1 at 



an + p 



(24) 



which points along the time evolution vector. Since we 
are considering orbiting binaries this Killing vector would 
have to be a helical Killing vector, which means that at 
spatial infinity (r — > oo) 



n a = T a and f3 a = 0$ Q 



(25) 



where T a and <& a are the asymptotic time-translation and 
rotational Killing vectors at spatial infinity and VL is the 
angular velocity with which the binary rotates. 
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A. Choice of lapse 

In order to allow for the helical Killing vector (|24() we 
now proceed by choosing a lapse such that 

d t K = (26) 

which leads to the elliptic equation 

V 2 (a<j>) - aV 2 <j) = o r ':n/x,,A" + (3 l V t K) (27) 

for the lapse. Note that this equation is valid for any 
and time independent K. We now use the specific form 
of puncture data and make the ansatz 



a<p - 



' c\m\ c 2 m 2 \ 
1 - I — h — I + v 



2n 



2r 2 



(28) 



for the lapse, where v is a finite correction. Together with 
Eqs. (|2*7ll and ljTT)l this ansatz yields the elliptic equation 
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V 2 w = -(a 



i LW ij LW ii , 



(29) 



which has to be solved numerically subject to the bound- 
ary condition 



lim v = 0. 

r — >oo 



Combining Eqs. (|2*U|) and yields 



V 2n 



2r 2 



in i 
2ri 



ni2 
2r 2 



(30) 



(31) 



we also get d t gij = d t <j> — 0. However d t Aij ^ 0, and 
hence c\ = — 1 is not a good choice since it leads to an 
evolving Aij even in the case of the static Schwarzschild 
spacetime. 

Another possibility obtained from c\ = 1 and v = is 



i _ m i 

1 _|_ Ull. ' 

1 T 2n 



(35) 



which also yields <9jA' = 0. This is the standard lapse of 
isotropic coordinates. If we again choose (3 l = we find 
that Eqs. i|23|l and (|24(l are now satisfied. Therefore c\ = 
1 is a good choice for a single BH. It turns out that for 
ci = 1 the Komar mass Ik (an, Soo) (see Eq. (131)1) below) 
computed from a equals the ADM mass of the BH. This 
is not true for c\ = — 1 for which Ik (an, Soo) = 0. 

Let us return to the general case of binary BH data and 
discuss what it means to impose equality of the Komar 
and ADM masses. The binary puncture data that we are 
considering may not possess a helical Killing vector 
but if they did then the Komar integral |2y, |23, 



^abcd 



(36) 



integrated over any closed 3-surface S containing the 
punctures would yield the same answer. Furthermore 
if we choose S to be a sphere at r — > oo and if £ a is 
normalized such that 



lim £ a n a = -1, 



(37) 



then 



Hence we obtain 



I K (tS) = M, 



ADM 



2fU 



ADM 



(38) 



lim a = 1 



(32) 



as boundary condition at spatial infinity, which is what 
one wants in asymptotically Minkowskian coordinates. 
The value at each puncture is 



is expected to hold as an immediate generalization of [2jJ • 
Here M^ DM and J^ DM are the ADM mass and angular 
momentum at spatial infinity (r — > oo) given by 



n.rADM _ pADM pADM 
oo ii OO V oo 



(39) 



lim a 

r A ^0 



~CA, 



(33) 



that is, our ansatz for the lapse l|28|l introduces the free- 
dom to pick the value of the lapse at the inner asymp- 
totically flat ends. 

This freedom turns out to be essential when trying to 
satisfy the mass condition for a helical Killing vector field. 
But let us first discuss the situation for a single puncture 
without momentum, i.e. the case of P[ = P| = ni2 = 
Kij = u = 0, which corresponds to a Schwarzschild BH 
in isotropic coordinates. Of course, the Schwarzschild 
spacetime is static and thus has a Killing vector £ a so 
that we can ask for which c\ Eqs. I)23|) and (|24|l do indeed 
hold. It is easy to see that 



a = 1 everywhere, 



(34) 



which corresponds to v — and c\ = —1, fulfills Eq. I|29|) 
so that dtK = holds. If in addition we choose j3 l = 0, 



and 



' D Qoo M — ^ J s (9ij,i 9ii,j)dS : ' 



P; 



ADM 



— \ (ICv - K VlJ )dSi 

o7T J c 

« <~>oo 



J 



ADM 



8vr 



(40) 



(41) 



(42) 



Note that the normalization condition (|37|l requires that 
a — ► 1 for r — > oo. 

If we insert £ a as given in Eq. (|24ll into Eq. I|36|l we 
obtain 

I K %S) = I K (an,S)+I K (f3,S) 

= TZ I ^i ad§i -JZ I KijP'dS 3 . (43) 



47T 



47T 
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Using Eq. I|25(l we find that the second term integrated 
at r — > oo is 



-2QJ r 



ADM 



(44) 

Therefore combining Eqs. (JSHJ, |@3J| and JSJ the condi- 
tion 



I K (cm ) S oc )=M£ DM (45) 

must hold if the helical Killing vector of Eq. I|24|) exists. 

We can now address the question what the values — ca 
of the lapse at the punctures should be. Since we are 
interested in the case where a helical Killing vector ex- 
ists (at least approximately), wc should pick ca such that 
Eq. 145JI is fulfilled, otherwise £° cannot be a Killing vec- 
tor. Since Eq. I|45|l is just one condition on in principle 
two unknowns, we also set c± = c 2 so that the lapse has 
the same value at each puncture. This is justified by the 
fact that in our numerical computations we only study 
equal mass binaries with m\ — m<i. Of course Eq. (|45|) is 
only a necessary condition for a Killing vector. To test 
whether a Killing vector really exists, we have to check 
how well Eq. (£3 is fulfilled. 

We summarize our numerical results in Tab.lTTl Given 
the ISCO and Pre-ISCO parameters described in Tab. [I] 
we iterate ca until I K (an 1 S 00 ) = M^ DM . Tab. lU gives 
the values for ca for which this can indeed be achieved. 
Since the lapse has now been fixed, we proceed to check 
two other relations, namely the mass equality at the 
punctures (as opposed to infinity) and the angular ve- 
locity predicted by the Komar integral. 

Note that Eq. I)38|) can be used to compute the angu- 
lar velocity fl if the Komar integral Ik(£, S) is evaluated 
on a surface S ^ Soo. Let us use Eq. (|43|l to calculate 
Ik(£,, S p ) for S p — Si U S2, where Si and S 2 are infinites- 
imally small spheres around each puncture. If we set (3 l 
to zero at the punctures (an assumption which will be 
justified below in Sec. IIII 5|l . the shift term in Eq. I)43|) 
does not contribute and we find 



kit, s p ) 



I K (an, Si) 

1 + Ci 



m 2 



2 

l + c 2 



I K (an,S 2 ) 

Vl + C\Ui 

2 

v 2 + C 2 U 2 



(ci - c 2 )m 2 



4L> 

(c 2 - ci)mi 



4L> 



(46) 

where ua and va are the values of u and v at puncture 
A. Note that ci = c 2 implies that there is no explicit 
dependence on D, but ua and va depend on D. If £ a is 
a Killing vector, Ik(S,, &p) should have the same value as 
in Eq. (|38|l . which in turn implies 



n = 



M™ M -Ik{S,S p ) 



(47) 



2J^m 

Comparing Tabs.lTlandlTTlwe see that the angular velocity 
(|47|l of the binary is very close to what is found with the 
effective potential method. 



TABLE II: The value of the lapse at the puncture for the 
three parameter sets of table [I] is chosen such that Eq. 1451 
holds. We find that then Eq. (1481 also approximately holds. 
In addition the £1 of Eq. 1471 is close to what the effective 
potential method predicts, and up to numerical accuracy Q = 
Q/3 for Q/3 introduced in Sec. 1111 51 



Parameter set 


ISCO 


Pre-ISCO 


Pa =0 


CA 


0.726 


0.829 


1.000 


M ADM 


1.003 


1.013 


1.000 


Ik (an, Sac) 


1.003 


1.013 


1.000 


Mi DM 


0.514 


0.518 


0.534 


Ii<(an, Sa) 


0.372 


0.428 


0.500 


c A Mi nM 


0.373 


0.430 


0.534 


c A M* um -I K {an,,S A ) 


0.0027 


0.0047 


0.0637 




M* DM Q. 


0.168 


0.100 


0.000 




0.168 


0.100 


0.000 



A further check is provided by the analogue 

I K (an,S A ) = c A M£ DM (48) 



to Eq. I|45|l . Here the ADM mass of each puncture is 
given by 



M ADM 



m A (l + u A ) + 



mim 2 
2L> : 



(49) 



and the factor ca comes from the fact that at the punc- 
ture £ a is normalized such that £ a n a = ca, which differs 
from the standard normalization in Eq. (|37|) . As one 
can see in Tab. ITflour numerical results for M^ DM and 
Ik (an, Sa) approximately satisfy Eq. (|48|l in the case of 
quasicircular orbits, while for the Pa — data Eq. fHfl) is 
violated. Yet even for quasicircular orbits the deviation 
from Eq. H48f) is larger than our numerical errors, but 
the parameters of Tab. [I] used to put the punctures into 
quasicircular orbits (according to the effective potential 
method) are not very accurate. We believe that these 
parameters are accurate up to 1% error, hence some de- 
viations from Eq. I|48|) are expected. In future work we 
plan to construct quasicircular orbits within our numeri- 
cal method, which involves varying P A , and it should be 
possible to define P\ for quasicircular orbits by imposing 
Eq. I|48|) . For the parameters used in the present work, 
we conclude that the one parameter freedom of the lapse 
at the puncture allows us to find a lapse that is approxi- 
mately compatible with the two necessary conditions l|45|) 
and H48fl for a helical Killing vector. 

A plot of the lapse is shown in Fig. |2 for two punc- 
tures located on the y-axis for the case of the Pre-ISCO 
parameter set of Tab.QJ 



Choice of shift 



In Sec. IIII Al we have chosen the lapse such that dtK = 
0. Similarly we would like to choose a shift which results 
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FIG. 2: The lapse a along the y-axis for the Pre-ISCO data 
set. The punctures are on the j/-axis at y = ±1.849. 

in dtgij — 0. In general (if the conformal factor is chosen 
such that dfdetgij = 0), we have 

dtgij = -2a<T 4 A 4 , + Lfoj. (50) 

Note however that L/3y is by definition purely longi- 
tudinal, so that the left hand side can only be zero if 
2a</>~ 4 Ay is purely longitudinal as well. The CTS con- 
struction enforces this feature by exchanging the role of 
Ay and dtgij , he. the time derivative of the metric is 
treated as free data that defines the extrinsic curvature 
through 

Aj = 2q *_ 4 {LPij - d t gij) . (51) 

For maximal slicing and vanishing time derivative of the 
metric one obtains an elliptic equation for the shift by 
taking a divergence of i|51[) and using the momentum con- 
straint which leads to 

VjL/3 ij = L(3 i: >V 3 log (ac/)- 6 ) . (52) 

In the CTS construction these are three of five coupled 
equations for 0, a, and /3\ 

However, for puncture data with maximal slicing, <f>, a, 
and Ay are determined independently of the shift, and 
in particular 2a</> Ay will not be purely longitudinal in 
general. In fact, for the puncture data we are considering 
here, 

2a<p- 4 A ij = 2a<j)~ <i LW l j 1 (53) 

so that 2ouj> Aij is the product of a longitudinal piece 
and a scalar function. 

But we can still determine a shift that completely re- 
moves the longitudinal piece from the time derivative of 



the metric. Concretely, in order to find the longitudinal 
piece Lfj^ of Eq. (|53|) for puncture data, let us make the 
ansatz 

2a<t>- Q LW lJ = Lfri + v^ T) (54) 

where vjl T is a possible transverse traceless piece with 

^ ' 3 v tt = 0- By taking the divergence of Eq. (|54|l we find 
the elliptic equation 

VjLF* = Vj (2ac/)- e LW lj ) , (55) 

which has a unique solution for (3 l for given boundary 
conditions (compare with Eq. lol'li for CTS). Once we 
have found /3 1 , we can then determine v^ T from Eq. 154fl . 
Since we are interested in a binary configuration in coro- 
tating coordinates we should adopt 

lim 0- = (56) 

r — >oo 

as boundary conditions for (3 l at spatial infinity. Here 
is the asymptotic rotational Killing vector pointing along 
the direction of rotation and ri^j is the angular velocity 
of the binary. Let us split the shift into the two pieces 

F=F + Pl ot , (57) 

where we have introduced the rotational piece 

0l ot = e>'fi* (58) 

and a piece with 

lim /3j = 0. (59) 

r — >oo 

The benefit of this split is that Lffi ot = 0, so that if 
Pq fulfills Eq. with boundary condition <|5§|) . f3 l of 
Eq. I|57|) immediately fulfills Eq. I|55|l with boundary con- 
dition (|56|l . We have numerically solved for /3q and then 
added a rotational piece of the form l|58|l such that (3 % — 
at each puncture. This is a very desirable property since 
punctures have no linear momentum when viewed from 
the asymptotically flat region at the puncture, so that 
the natural choice for the shift is indeed zero. The shift 
for the Pre-ISCO parameter set of Tab. [I] is shown in 
Fig. |3| Notice that f3 l is linear for large r in agreement 
with boundary condition l|5()|) . It turns out that up to 
the accuracy of our numerical computation, the re- 
quired to achieve (3 % — at each puncture is equal to the 
O computed from Eq. I|47l) . We want to stress here that 
fig = O is another necessary condition for a Killing vec- 
tor (see appendix^) and not expected to hold a priori 
for puncture data. 

In Fig. 0| we show plots of the three terms in Eq. H54fl . 
We find that vjl T in Eq. i|54|) is indeed non-zero. Hence 
Eq. I|5U|) becomes 

dtgij = -2a<t>' Q LW l3 + Lfcj = -v^ T , (60) 
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FIG. 3: The shift f3 x for Pre-ISCO data along the y-axis, 
which connects the two punctures. At each puncture /3 X van- 
ishes. 
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FIG. 4: The three terms in Eq. E3) for Pre-ISCO data. 



and dtgij = cannot be achieved for puncture data. 
Note however that v^ T is smaller than , and that 
2a<t>" & LW 1 ^ is purely longitudinal for large r. Thus dtgij 
is at least reduced by the above choice of f3\ The maxi- 
mum value is 



ADM 



flu 



(61) 
evolves is 



max d t gij = 2x 10~ 2 /M 

Hence the smallest timescale on which 
T 9ij = 27r(l/max<9 t g y ) - 300M ADM which is longer 
than the dynamical timescale given by T^ yn = 2tt/Q, ~ 
60 M ADM . This means that is approximately con- 
stant up to an error of at most Td yn /T gij ~ 20% in the 



center and much less away from the punctures. 

So far we have only addressed the time derivatives of K 
and gij. We will now use the BSSN evolution equations 
to investigate the remaining time derivatives. Recall that 
BSSN introduce the additional variables 



I ' = <)/ ■>. 



(62) 



8ij on the initial slice 



where the BSSN metric 7^- = g^ 
for the case of conformally flat puncture data. Using 
Eqs. (J5DJ and we obtain 



d t r = -Vjd t g ij = V j Lf3 ij - Vj (2a(f>- 6 LW lj ) . (63) 

Thus, when Eq. I|55|) is fulfilled we have 

dtf 1 = 0, (64) 

so that Y % does not evolve with our choice of shift. 
Put differently we have three free functions available in 
terms of the shift vector. This does not suffice to set 
d t gij = for puncture data, but we can impose the so- 
called Gamma- freezing condition l|64[l [rH |23| , which we 
have obtained here by constructing a shift for CTT data 
motivated by CTS data. 

The time derivatives of the remaining variables have 
to be studied numerically. We will do this next for the 
Pre-ISCO parameter set of Tab. |U First we look at the 
variable 



ip = log < 



(65) 



which BSSN introduce by decomposing the physical met- 
ric as gij — e ffyij. Numerically we find that dtf is very 
close to zero. The maximum value is 



max d t f = max 



d t 4> 



5 x 10" 



(66) 



This means that <fi evolves on a timescale T p hi — 
2ir((f>/d t (t>) > 1200M ADM , which is much longer than the 
dynamical timescale given by Td yn ~ 60M ADM . 

Finally, Fig. shows the BSSN variable Ay = 



Aij, and its time derivative. We see that 



dtAij does not vanish. In fact for some components such 
as A yy , dtA^ is not much smaller then Aij itself, while 
the time derivative of the dominant component A xy is 
small. This implies that puncture data cannot have a 
helical Killing vector. However as we can see in Fig. El 
dtA^ is reduced by our choice of lapse and shift when 
compared to the case of a = 1 and (3 l — 0, which is of- 
ten used as initial data for lapse and shift when evolving 
with dynamical gauge conditions. Hence our choice of 
lapse and shift at least brings us closer to the situation 
of minimal evolution. 



IV. DISCUSSION 

We have shown that puncture data for quasicircular 
binary BHs obtained in the CTT framework with the 
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FIG. 5: Selected components of Aij and their time derivatives 
along the diagonal x = y — z for Pre-ISCO data. 
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FIG. 6: The time derivative of A yy for a and /3 ! computed 
from Eqs. 1271 and 1551 is much smaller than for a — 1 and 
ft = for Pre-ISCO data. 



effective potential method allows a special gauge choice 
that realizes some of the necessary conditions for the ex- 
istence of an approximate helical Killing vector field. In- 
troducing a free parameter for the lapse at the punctures 
we can satisfy the condition that the Komar and ADM 
mass agree at spatial infinity. Since there are no fur- 
ther free parameters in our gauge choice, it becomes a 
non-trivial fact that the Komar and ADM mass also ap- 
proximately agree at the punctures, and that the angular 
velocity given by three different methods agree within a 



certain accuracy. The angular velocity ftp determined by 
the shift that vanishes at the punctures is equal to the 
angular velocity fl inferred by assuming that the Komar 
integral can be evaluated on any surface containing the 
punctures, and both Qp and agree approximately with 
the angular velocity prescribed by the effective potential 
method. 

Our gauge choice controls some of the time deriva- 
tives of the metric variables, dtK — and dtdjg 1 ^ = 0. 
Furthermore, we find that <f> varies on a timescale much 
larger than the orbital timescale and that gij varies on a 
timescale which is at least five times the orbital timescale. 
The magnitude of dtAij is reduced when we compare to 
the case of lapse a = 1 and zero shift, but it still varies 
on the orbital timescale. 

One conclusion to draw from our results is that qua- 
sicircular puncture data is not as far from possessing a 
helical Killing vector as one might have expected. This is 
an important observation since to date all gravitational 
wave forms obtained for binary BH mergers are based on 
such initial data. 

Note that in the present work we have used parame- 
ters for quasicircular puncture data found in the litera- 
ture. It should now be possible to construct sequences 
of quasicircular orbits based on the Komar and ADM 
mass equalities and compare with the effective potential 
method. At a given separation D, one has to vary the 
linear momenta P\ and the constants ca for the lapse at 
the punctures until both Eqs. (|45|l and i|48[) are satisfied. 
The variation to find quasicircular orbits has to keep, 
for example, the ADM masses at the punctures constant 
(cmp. 30]), which would be particularly natural in our 
approach. 

Another direction for future work is to perform evolu- 
tions of puncture data with the initial gauge derived here. 
Note that this is not entirely straight-forward because the 
lapse crosses from positive to negative values near the ap- 
parent horizon, which is also true for the Meudon data 
sets construct in the CTS framework, but in that case 
an isometry condition is imposed. It is not clear whether 
there are numerical problems related to a negative lapse, 
even if one applies BH excision techniques, but this is 
certainly worth exploring. 

Finally, an immediate question is whether one can per- 
haps modify the extrinsic curvature of the puncture data 
to better control <9t Ay . Puncture data provides a genuine 
technical simplification over methods with excision, but 
so far it has not been possible to construct thin sandwich 
puncture data in order to improve the extrinsic curva- 
ture. One aspect of the problem becomes already appar- 
ent when considering maximal slicing of Schwarzschild, 
which is known analytically 31] . It is not hard to see 
that demanding dt§ij = for all t immediately leads to 
the standard isotropic coordinates with negative lapse 
at the puncture and vanishing shift 32]. The conclu- 
sion is that there is no maximal slicing of Schwarzschild 
with static metric components and positive lapse. On 
the other hand, in evolutions of Schwarzschild with 
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positive lapse have been performed that display an ap- 
proximate Killing vector, therefore the existence of thin 
sandwich puncture data has not been ruled out. 
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1 

4-7T 



6d t log 



A»rd t9ij - d t K 
2 



(A5) 



Next, let us define Ik (an + (3, S) by the right hand side 
of Eq. l|43|l even if no Killing vector £ exists. Now, if the 
volume V in Eq. I|A5J| is bounded by the surfaces 6*00 and 



S p we find that 



I K (cm + (3, Soo) - I K {an + (3, S p ) = ly, 



(A6) 



APPENDIX A: INTEGRAL I K DEFINED 
WITHOUT A KILLING VECTOR 

If the conformal factor is chosen such that dtdetgij = 0, 
the evolution equations for </>, gij and K are 

6d t log0 = -aK + Vt(3\ (Al) 
dtgij = 0- 4 (L/3 l3 -2aA lJ ) (A2) 
d t K = V 2 a + a4A !J +aK 2 + (3 l d l K. (A3) 

If we solve Eqs. I|A1I) and l)A2|l for aK and A^- respec- 
tively and insert the results into Eq. I|A3|I we obtain 

^(Via-Kij^) = -A ij <t>% 9ij -d t K 



65 t log^--V l /3 l ) K,(AA) 



where we have made use of the momentum constraint. 
Now note that the term in parentheses on the left hand 
side is the same as in Eq. (|43H . This suggests to define 
the volume integral 



1 



i.e. the value of Ik{coi + P,S) in general does depend 
on the surface S on which it is evaluated. If, however, 
an a +P a is a Killing vector, ly vanishes and /R-(an+/3, S) 
does not depend on S. 

We now adjust the shift such that (3 l — on S p and 
(¥ = 51/3$' on Soc, then Eq. (|A6|) becomes 

I K (an, Soo) - 2Hp J* DM - I K (an, S p ) = ly (A7) 
Choosing the lapse such that Eq. (|45|1 is satisfied, yields 



M. 



ADM 



Ik{£, S p ) - 1\ 



2J£P M 



(A8) 



Hence Hp and of Eq. I|47() differ by the term ly, which 
vanishes if an a + f3 a is a Killing vector. Another neces- 
sary condition for the existence of a Killing vector is thus 
Qp = f2. For punctures with our choice of lapse we have 
K = dtK = 0, so that all terms involving K in ly of 
Eq. (|A5|1 vanish, the term with dtgij in Iv however may 
not integrate to zero, so that Hp = Q is non-trivial for 
puncture data. 
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